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We investigate the origin of the depolarization rates in ultrathin adsorbate-stabilized ferroelec¬ 
tric wires. By applying density functional theory calculations and analytic modeling, we demon¬ 
strate that the depolarization results from the leakage of charges stored at the surface adsorbates, 
which play an important role in the polarization stabilization. The depolarization speed varies 
with thickness and temperature, following several complex trends. A comprehensive physical model 
is presented, in which quantum tunneling, Schottky emission and temperature dependent electron 
mobility are taken into consideration. This model simulates experimental results, validating the 
physical mechanism. We also expect that this improved tunneling-Schottky emission model could 
be applied to predict the retention time of polarization and the leakage current for various ferro¬ 
electric materials with different thicknesses and temperatures. 


I. INTRODUCTION 

Spontaneous electric polarization makes perovskite- 
based oxides of great interest for application to non¬ 
volatile memory devices [1, :2|. However, the polarization 
of ferroelectric materials may not be infinitely stable, and 
retention time is one of the key factors determining the 
performance of memory devices in nonvolatile technol¬ 
ogy. The proposed reasons for the polarization instabil¬ 
ity have included the depolarization field and the leak¬ 
age current [d|, whose effects become more significant 
as the oxide film gets thinner. Therefore, for success¬ 
ful technology application, the depolarization processes 
of nanoscale ferroelectric oxides must be better under¬ 
stood. Here, we report a combined experimental and 
theoretical investigation of the depolarization process of 
single-crystalline BaTiOs nanowires. We attribute the 
decay of polarization to the leakage of surface screening 
charge and propose an analytical model to explain the 
experimental decay rates. 


A. Experimental background 

The effects of the depolarization field on the stability 
of ferroelectricity in ultrathin materials was explored in 
BaTiOs nanowires by measuring the ferroelectric tran¬ 
sition temperature as a function of the nanowire diam¬ 
eter in the range of 3-48 nrn &E3- Positive ferroelec¬ 


tric domains were written perpendicular to the nanowire 
axis using a negative-bias voltage (—10 V) applied by 
a conductive scanning probe microscope cantilever tip 
(under ultrahigh vacuum conditions with a base pres¬ 
sure of 10“ 10 torr). The time evolution of the polarized 
domain was then monitored via time-resolved measure¬ 
ments of the local electric field, using non-contact elec¬ 
trostatic force microscopy (EFM). The writing and read¬ 
ing processes were done at various temperatures, starting 
from ps 393 K for thin nanowires (3-11 mn) and ps 418 
K for the thicker ones (12-37 mn). The two sets of 
nanowires were progressively cooled down and retested 
to ps 308 K and « 383 K were reached respectively. The 
Curie temperature Tq was defined as the highest tem¬ 
perature below which the polarization signal persists for 
a period longer than 200 hours. Experiments showed 
that Tc is inversely proportional to the diameter of the 
BaTiC >3 nanowire, in accord with standard models of de¬ 
polarization field [8j. At several temperatures above Tc, 
the surface potential signals were measured during the 
process of polarization decay. The magnitude of surface 
potential was fitted with the expression: 

S(t) = S( 0)e~ kdt (1) 

where S(t ) is the potential at time t, which is propor¬ 
tional to the surface screening charge, and kd is the decay 
rate, with the unit s _1 [§]. Here, we should note that ex¬ 
perimentally it is found that the signals decay with time 
approximately (but not perfectly) exponentially. Despite 
the slight deviations, the decay rate kd obtained from the 
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TABLE I. Decay rates of ferroelectric polarization for dif¬ 
ferent BaTiC >3 nanowire thicknesses and temperatures (thin 
nanowires). 


Temperature 

Decay Rates kd (s 1 ) 

5 nm 

7 nm 

9 nm 

393.6 K 

0.0287 

0.00994 

0.00368 

389.4 I< 

* 

0.0084 

0.0027 

384.3 K 

0.0232 

0.00485 

0.00162 

379.5 I< 

0.0161 

0.00248 

0.000884 

373.9 K 

0.00507 

0.00132 

0.0001 

369.3 K 

0.00301 

0.000413 

* 

364.8 I< 

0.00192 

* 

* 

359.7 K 

0.000916 

* 

* 

354.5 I< 

0.000287 

* 

* 

349.2 K 

0.000267 

* 

* 


* Data not available 


TABLE II. Decay rates of ferroelectric polarization for dif¬ 
ferent BaTiC >3 nanowire thicknesses and temperatures (thick 
nanowires). 


Temperature 

Decay Rates kd (s x ) 

21 nm 

25 nm 

37 nm 

48 nm 

419 K 

0.00333 

0.00263 

* 

0.00171 

410 I< 

* 

0.00122 

0.0021 

0.00197 

405 K 

0.000989 

0.000997 

0.00084 

0.00105 

399 I< 

0.000244 

0.00027 

* 

* 

395 K 

* 

0.000114 

0.000115 

* 


* Data not available 


data fitting is still an important physical parameter in 
describing the polarization decay speed. From the time 
evolution of the signal [5j], we see that the bright circu¬ 
lar signal faded without expansion, which means that the 
depolarization is a process of leakage or tunneling, rather 
than diffusion, of the surface screening charge. The ob¬ 
served decay rates kd , which vary with nanowire thickness 
and temperature, are presented in Tables I and II. 

The experimental data show three general trends: (1) 
The depolarization process is slow (several hours); (2) 
For any thickness, the decay rate kd increases with tem¬ 
perature; (3) For thin nanowires (5-9 nm), the decay 
rate kd changes with thickness dramatically. However, 
for thicker wires, kd stays nearly constant at different 
thicknesses (21, 25, 37 and 48 nm). Our study aims at 
illustrating the physical essence of these trends. 

B. Polarization stabilization 

In recent years, many studies have investigated the de¬ 
pendence of polarization stabilization and leakage current 
on chemical environment, temperature, electrode mate¬ 
rial, and thickness [I |, HolSj . From these studies, 
several basic principles could be drawn: 

(1) Surface polarization charge should be compensated 
by screening charge, in order to passivate the depo¬ 


larization field and stabilize the ferroelectric distor¬ 
tion, or else the polarization would become unsta¬ 
ble. The screening charge could be stored in surface 
electrodes or adsorbates @, ITol - fT8j ; 

(2) The polarization state of the material may lead to 
preferential adsorption of certain molecules on the 
surface [Hj-[25|; 

(3) The response of polarization with electric field or 
temperature is fast, but the time scale for dissipa¬ 
tion of the surface screening charge is slow (hours 
or days) (2(14281 . 

Based on the evidence above and the observations in 
our experiments, we propose that the physical process 
of depolarization in the BaTiCU nanowire experiment is 
as follows: after the polarization is written, surface ad¬ 
sorbates on the nanowire act as an electrode that stores 
screening charge and stabilizes the polarization p|. For 
the case without external applied voltage and above Tc, 
the polarized state is not stable. But due to the sta¬ 
bilization of screening charge in the surface adsorbates, 
polarization in the nanowire still persists for some time. 
Screening charge leaks from the top electrode (surface ad¬ 
sorbate) to the other side of the BaTiCU nanowire (gold 
substrate). At the same time, polarization reduces along 
with the screening charge. This process is slow and takes 
hours. 

In the following parts of this paper, support, analysis, 
and modeling of the physical processes described above 
are shown. In Section II, we use density functional theory 
to demonstrate the role of surface adsorption in surface 
charge screening. In Section III, analytical expressions 
describing leakage current leading to depolarization are 
developed. Finally, in Section IV, we present the results 
and discussion. 


II. DENSITY FUNCTIONAL THEORY 
CALCULATION 

In order to construct a theory of the depolarization 
process, density functional theory calculations (DFT) are 
carried out to assess the role of surface molecular and 
atomic adsorbates. We investigate the OH molecule on 
BaO-terminated BaTiCU slabs, as OH is the predominant 
species found on oxide surfaces, as demonstrated by both 
infrared spectra and ab initio calculations pi. [29I.I30]. 

A supercell slab method was used to describe the sys¬ 
tem. The supercell consists of the BaTiC >3 slab with 
11 atomic layers, a full coverage of OH (one adsorbate 
per adsorption site) on the lxl BaTiOs surface, and a 
vacuum of more than 20 A. The atoms are represented 
by norm-conserving pseudopotentials generated using the 
OPIUM [3l| code with a plane-wave cutoff of 50 llvlv4 . 

From the relaxed structure, we see that the presence of 
the OH adsorbates enhances ferroelectricity at the pos¬ 
itively polarized surface and maintains a characteristic 
ferroelectric displacement pattern throughout the film. 
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FIG. 1. Relaxed structures for the BaTiOs/OH systems ac¬ 
quired from DFT calculation. Green, blue, red and gray 
spheres represent barium, titanium, oxygen and hydrogen 
atoms respectively, (a) Polar system, OH is adsorbed at the 
surface; (b) Non-polar structure, OH is far away from the sur¬ 
face of the BaTiOa surface; (c) Top view of the OH adsorbate 
and the first atomic layer of BaTiOa. 


TABLE III. Lowdin population in the orbitals of OH adsorbed 
to BaO-terminated BaTiOs. 



Polar 

Non-polar 

H Is orbital 

0.6121 

0.5601 

0 2s orbital 

1.7315 

1.8138 

O 2 p orbital 

5.0719 

4.5980 

Net Charge 

-0.4155 

0.0289 


The density of states projected onto atomic orbitals 
(PDOS) was calculated to characterize the charge distri¬ 
bution on each atom. Results are shown in Table III, 
from which we see that if the BaTiOa nanowire is posi¬ 
tively polarized, the hydroxyl oxygen 2 p orbital possesses 
more electrons and OH is overall negatively charged. 
This is a robust evidence demonstrating that surface ad¬ 
sorbates stabilize the polarization by holding screening 
charges, which has an effect similar to an electrode. The 
leakage of charge mainly from the 2 p orbital of oxygen 
then results in the decay of polarization. 


III. FORMALISM OF LEAKAGE CURRENT 
CALCULATION 

There have been previous reports calculating the mag¬ 
nitude of leakage through ferroelectric films with tunnel¬ 
ing models, Schottky emission, and the modified Schot- 
tky equation |33l - 135j . However, all these models come 
across difficulties in explaining all three depolarization 
rate trends mentioned in Section I. B. A tunneling model 
alone cannot explain that for thick nanowires, decay rates 
are nearly thickness independent and for any thickness, 
the decay rates are strongly temperature dependent. On 
the other hand, a Schottky emission model cannot ac¬ 
count for the thin nanowire thickness dependent rates. 
Here, we build a comprehensive model from the “effec¬ 
tive velocity” point of view, which both accounts for the 
experimental results and illustrates their physical mech¬ 
anisms. 

The general expression for the time ( t ) dependent leak¬ 
age current J can be written as 

~ = J (t) = J Q (t)n (fc) v e s (fc) d 3 fc (2) 

n (fc) is the probability that an electron possesses a wave 
vector between k and k + dk. v e g ( k ) is the effective 
velocity along the z direction, which is normal to the 
surface of the BaTiOa nanowire, for the electrons with 
wave vector k. Q is the amount of extra charge (com¬ 
pared with neutral OH) stored in the adsorbate. In the 
following subsections, we define the parameterization of 
equation ©. 


A. Wave vector distribution of electrons in 
adsorbate 


Unlike in traditional metal electrodes, electrons occu¬ 
pying orbitals localized on the OH adsorbates cannot be 
treated as a free electron gas, and the wave vector dis¬ 
tribution does not follow Fermi-Dirac statistics. Instead, 
the wave vector spectrum can be estimated from Bessel- 
Fourier transformation of the 2 p orbital of oxygen, since 
screening charge is mainly associated with this atomic 
orbital. Here, the radial part of the 2 p orbital of the 
oxygen is represented by a double-zeta function [3f|. 


4>2 P (r) = R (r) Yjo (6 r , <j> r ) 


= c il 


( 2 *i r 


4! 


re~ zir + c 2 1 


'(2 Z2 )" 


4! 


re~ Z2r Yiq (0 r ,0 r ) 


= (c' 1 re~ Zir + c' 2 re~ Z2r ) Y w (6 r , tj> r ) 

<Mfe) = - 4/2 [ e~ lk r K(r)d 3 

(27T) J 


— — \j —4tFi 0 (0fc, (j>k) 


d x z\k c' 2 Z 2 k 


(zf + fc 2 ) (z 2 + fc 2 ) 


( 3 ) 

( 4 ) 
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n{k) = \(j) 2p (fc)| 

1 98 

= — 1^10 (d k Ak)\ 2 


c^zik 


c' 2 Z2k 


(zf + k 2 ) (z 2 + k 2 ) 


( 5 ) 


Yio is the spherical harmonic for l = 1 and m = 0 . 
ci _2 and zi t 2 are the parameters in the double-zeta func¬ 
tion, acquired from previous reference [36!]. and c' 2 
are reduced coefficients taking the normalization factors 

\J *• 2 ] 2\ and \J ^ 4 ^ into consideration. 0 Ty k and <fi rt k are 
the angles between the directions of r, k and the axes in 
spherical coordinates. </> 2p (r) and (j) 2 p ( k ) are expressions 
for the oxygen 2p wavefunction in coordinate and wave 
vector representations. In this way, we obtain an analyti¬ 
cal expression for n{k). Here, we should note that in part 
II, we used pseudo wavefunctions in the DFT calcula¬ 
tions; compared with all-electron wavefunctions, pseudo 
wavefunctions have lower high-fc components on purpose 
to limit the number of plane waves used [3?j . This does 
not affect the accuracy of charge leakage rate calculation. 
This is because the DFT calculation is used only to illus¬ 
trate the role of screening charge and the mechanism of 
the depolarization process. The wave-vector distribution 
in this model is derived from the double-zeta function as 
described above. Additionally, in the later discussion, we 
will also include the fact that high-speed electrons lose 
their initial momentum quickly and drift under the ef¬ 
fect of the electric field. Therefore, an underestimation 
of the high-fc components has little effect on the charge 
dissipation speed. 


the wave vector pointing toward the nanowire k ■ z > 0 . 
Otherwise, the electron does not contribute to the leak¬ 
age current J and its effective velocity u e fi = 0. For the 
electron moving toward the nanowire, the expression for 
effective velocity varies depending on whether the energy 
of the electron is higher than that of BaTiOs conduction 
band edge. The potential energy of an electron affiliated 
with the hydroxyl molecule oxygen 2p orbital is set as Vo . 
The total energy for this electron can be written as the 
sum of its kinetic energy T {k) and Vo: 

E(k) = T(k) + V 0 = ^- + V 0 . ( 6 ) 

2 m o 

As shown in FIG. 2, the conduction band for the BaTiC >3 
nanowire is not flat, and the slope equals the electronic 
charge e times the electric field Ef e inside the ferroelectric 
nanowire. Therefore, 


V 2 = V 1 - eE te d 0 , (7) 

For the case E (k) < V 2 , the mechanism that governs the 
electron movement is quantum tunneling. The electron 
tunnels from the surface adsorbate hydroxyl to the gold 
electrode at the other side, through the BaTiC >3 nanowire 
as an energy barrier. Around Tq , the dielectric constant 
of BaTiC >3 is large and E{ e is small. Therefore, the con¬ 
duction band is nearly flat and we use the approximation 
that the energy barrier is a cuboid with the height V 2 . 
The transmission coefficient P (fc), which is also the prob¬ 
ability of one electron with wave vector k penetrating the 
barrier, could be expressed [38]: 


B. Effective velocity 

Here, we present an effective velocity model of 
the charge dissipation. The band diagram of the 
OH/BaTiOs/Au substrate system is shown in FIG.2. We 



V 0 


P(k) = 

4 


4 + 


[m 0 E + m* (V 2 - E )] 2 . 2 

mom*E (V 2 — E) 


' 2 m*d 2 (H 2 -^)] 1/2 

V 2 

( 8 ) 


where d is the length of the barrier, m* is the effective 
mass of electrons in BaTiOs- Assuming that the angle 
between the incident direction of the electron and the 
normal direction of the BaTiOs nanowire is 9 , as shown 
in FIG. 3, d could be expressed as 


do 

cos 9 


(9) 


In this tunneling process, an electron with an initial 
velocity v — ^ has a probability P (k) of passing the 
energy barrier with the thickness d. Thus, the effective 
velocity could be expressed as 


OH BaTi0 3 Au 

FIG. 2. Band diagram of the Au electrode/oxide insula¬ 
tor/adsorbate electrode system. 

only consider the effective velocities of the electrons with 


v eS (k) = vP (k) = * P(k) = fikc0S ^ P(k) (10) 

mo mo 

When the hydroxyl electron energy is higher than that 
of the BaTi0 3 conduction band edge E (k) > Vi, the 
electron is not classically forbidden to enter the BaTiOs 
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FIG. 3. Relationship of barrier thickness and direction of 
wave vector. 


nanowire. However, very high velocity electrons could 
have very short mean free paths. Moreover, in the 
BaTiC >3 crystal, which is an insulator and possesses a 
small mean free path (39|, the electron loses its initial 
momentum quickly and drifts under the effect of electric 
held £fe inside the ferroelectric nanowire: 

"Cell — I^Ef e , (H) 

where /r is the electronic mobility, which depends on 
temperature, intermediate energy levels and lattice vi¬ 
brations. Lattice vibrations have a significant influence 
on the drift mobility due to electron-phonon interac¬ 
tions [40l - [43l |. in which an electron may lose or gain en¬ 
ergy in a collision. This effect will be discussed in the 
following subsection. 

For the final case, V 2 < E ( k ) < V\, the electron first 
tunnels through an energy barrier and then drifts in the 
conduction band. The effective velocity is given by 

v eS = nE (e P (k). (12) 

For this case, the energy barrier goes to zero within the 
nanowire. At this the condition, the Wentzel-Kramers- 
Brillouin (WKB) [dij approximation fails. Since E (k) is 
close to Vi, we approximate P (k) to 1. 

Even though the approximations in the first and last 
cases overestimate the tunneling current, later we will 
point out that despite this, the tunneling current is still 
a negligible contribution to the overall current. 


C. Electric Properties of BaTiCE Nanowire 


The relationship of surface charge density Q, electric 
displacement D[ e and electric field Ef e th oug h the insu¬ 
lator is given by dielectric constant e ( T ) jlH ]: 


Q ( t ) = D fe = e (T) e 0 E te 


(13) 


According to the Lyddane-Sachs-Teller |45| relation, 
the temperature dependent static dielectric constant is 
expressed as: 


e(T) n ^t(T) 

£00 EL w Ti ( t ) 


(14) 


w is the frequency of lattice vibration and T,; and L,; rep¬ 
resent transverse and longitudinal modes. Approaching 
Tc , a transverse optical mode wto becomes “soft”, indi¬ 
cating the occurrence of phase transition [46l - l48j : 

u 2 to ^A{T-T c ) (15) 


cpy 1 cx^ 0 = A'(T-T c ) (16) 

A and A' are constants, and for BaTiC >3 bulk jH, 

Tc ~ 393 K. However, for a thin film or nanowire, if 
polarization exists, a depolarization field is induced due 
to the incomplete charge compensation of polarization 
charge [H, HH, HU- The depolarization field, which is 
anti-parallel to polarization, becomes significant as the 
thickness decreases. The depolarization field applies an 
electric force on the ions, and as a consequence, the 
soft vibrational mode is hardened. With the depolar¬ 
ization field effect included, the temperature dependent 
soft mode frequency for thin films should be rewritten as 


ft: A(T — T'c) 

(17) 

e (T) -1 = A' (T — T' c ) 

(18) 

T'c <T C 

(19) 


The depolarization field effect hardens the soft mode, in¬ 
creases the vibrational frequency, suppresses the ferro- 
electricity and lowers the Tc- The values of T' c for dif¬ 
ferent thicknesses could be inferred from Tc, which were 
previously experimentally measured Q . 

Previous studies illustrated that the electronic mobility 
in BaTiC >3 varies over a large range (10 -3 -10 _1 cm 2 /V- 
s) [H, Hif l54j . The variation of mobility could be at¬ 
tributed to the different trap levels, types and concen¬ 
trations from the fabrication process. Despite the uncer¬ 
tainty, there is a general rule that electron mobility is 
determined by the optical phonons [iol - |43j . 

Many studies in recent years pointed out that scat¬ 
tering by the soft transverse optical mode is the pri mary 
factor affecting electron drift ^j,[55|,[5a] . In this case |5g]: 

M = oc/(T)w ™ (20) 

and / (T) depends weakly on T for BaTiC >3 f55j,[56}. For 
T>T \£, 

b « Bu4o =B'(T- T'c ) (21) 

B and B' are constants. An explicit explanation about 
the soft mode dependent mobility was proposed in 
Ref. [43|. In brief, according to thermodynamics, the 
mean squared polarization fluctuation SP is related to 
the dielectric constant as [57]. 

(SP 2 ) = k B Te/V (x l/o-4 0 


( 22 ) 
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ks is Boltzmann constant and V is the volume. For the 
case T > Tq, a high soft mode frequency means a small 
dielectric constant, large mean-square TO phonon am¬ 
plitude and polarization fluctuation. Besides, in a small 
region with an approximately uniform polarization fluc¬ 
tuations, shifts in the conduction band edge A£ c is [58] 

A£ c = const + (35P 2 (23) 

/? is the polarization-potential parameter, which has a 
typical value j3 ~ 2 eV m 4 /C 2 . In the nanowire, the 
low soft mode frequency leads to inhomogeneity in A£ C1 
which then results to locally different effective masses 
and electronic energies. Therefore, an electron traveling 
through the wire scatters more. This is similar to an elec¬ 
tron traveling on a curved path and harder to accelerate. 
A lower soft mode frequency means more scattering, a 
shorter relaxation time, r e , and a smaller electron mobil¬ 
ity, since electron mobility is given by [59J 


Thus, in this simulation, we use the empirical expression 
of electron mobility as shown in equation m- 


of the nanowire. 

Vi - E(k) « 3.20 eV (26) 

where E (k) is the average energy of electrons in the ad¬ 
sorbate oxygen 2 p orbital, and 

' Ejk) = Tjkj + V 0 (27) 

_ r°° tffo 2 

T (k) = / n(k)~ - d 3 k = 68.28 eV (28) 

Jo 2 m 0 


Vl -Vo = Vl -E(k) + T(k) = 71.48 eV (29) 

In the modeling, the surface charge density Q decays 
with time t, but not exactly exponentially. We fit each 
Q vs. t curve with an exponential function by the least 
squares fitting method. In this way, we obtain decay con¬ 
stants from this model that can be compared with exper¬ 
imental ones. B' of equation ED is the only parameter 
calculated by fitting the data in experiments to the de¬ 
cay rates calculated by this analytical model, rather than 
from any references or ab initio calculation. 


IV. RESULTS AND DISCUSSIONS 


TABLE IV. Parameters involved in the presented tunneling 
and modified Schottky model of nanowire depolarization. 


With the expressions of effective velocity deduced in 
different wave vector range, 


if E (k) < V 2 , k ■ z > 0 
if V 2 < E(k) < Vi, k ■£ > 0 
if E(k) > Vi, k ■» > 0 
if k ■ z < 0, 

(25) 

we could calculate the time evolution of surface charge 
with equation ©. 

The expression of current density was shown in equa¬ 
tion (1). In the calculation of charge density evolution 
with time, the initial charge density used in the simula¬ 
tion is Q(0) = 0.26 C/m 2 [33], which is the spontaneous 
polarization for the tetragonal phase of BaTiC> 3 . The 
time window is selected as 10 4 s, which is long enough to 
demonstrate the general trend of time evolution of sur¬ 
face charge in adsorbate OH. Other parameters involved 
in the presented tunneling and modified Schottky model 
are listed in TABLE IV. Most of the parameters are 
from previous references. The difference of the potential 
energy of electrons in the adsorbate oxygen 2 p orbital 
(Vo) and in the conduction band formed by titanium 3d 
orbitals (Id) is estimated from the band gap of BaTiC >3 
(3.20 eV) [601. This is a good approximation because the 
adsorbate oxygen 2 p orbital is approximately at the same 
level with the valence band formed by oxygen 2 p orbitals 


UefT (k) = < 


— ■P (k) 

nE ie P (k) 
P Efe 
0 


Parameter Description 


T'c (5 nm) 

Tc for 5 nm nanowire a 

T'c (7 nm) 

Tc for 7 nm nanowire“ 

T'c (9 nm) 

Tc for 9 nm nanowire a 

T c 

Tc for thick nanowires a 

Cl 

Parameter in equation ([3]) ! 

Cl 

Parameter in equation (pH) { 

Zl 

Parameter in equation (|3|) ! 

Z2 

Parameter in equation (|3|) ! 

m* 

Effective electronic mass 

Vi-Vo 

Energy barrier d 

A' 

Defined in equation (1181) e 

B' 

Defined in equation (1211) ^ 


Value 
340.3 I< 

355.2 IC 
367.9 I< 

391 I< 

0.72540 

0.35173 

1.62807 bohr -1 
3.57388 bohr -1 
6.5 mo 
71.48 eV 
7.84X10 -6 K -1 
5.025xl0 -6 


“ Reference 4| 

6 Reference 36] 
c Reference 53| 

d From the band gap estimation 
e Reference [13] 
f The unit is cm 2 /VsK 


The decay rates calculated from this model and exper¬ 
imental data are shown in FIG. 4. Experimental data 
are marked with red circles and theoretically calculated 
ones are blue squares. Symbols for decay rates of a given 
thickness are connected with dashed lines, to demon¬ 
strate their temperature dependence. From the compari¬ 
son, it could been seen that this model not only simulates 
the experimental results to a good extent, but also sheds 
light on the general rules of ferroelectric leakage current. 

For nanowires with any thicknesses, calculations based 
on equations © and ED demonstrate that tunneling 
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Depolarization Rates 

Thin nanowires 


o- -o Experimental 
o- -n Theoretical 




FIG. 4. Comparision of experimental and theoretical depo¬ 
larization rates calculated from the presented model. 


current provides a negligible part of the total leakage 
current. Even for the thinnest nanowires (5 nm) with a 
surface charge density Q (0) = 0.26 C/m 2 , the tunneling 
current is around 10~ 6 A/m 2 or less, generally no more 
than 1 percent of the total leakage current. Emission 
of electrons with energies higher than the barrier con¬ 
tributes nearly all the current. The fact that tunneling, 
whose rate is thickness dependent, plays a minor role is 
consistent with the experimental finding that decay rates 
are thickness independent for thick nanowires & 

For thin nanowires, at a certain temperature which is 
above Tq, the soft mode frequency decreases with thick¬ 
ness due to the depolarization field effect. A thinner 
nanowire corresponds to a higher soft mode frequency 
and a smaller dielectric constant, as shown in equa¬ 
tion m- Thus, the electric field through a thinner 
nanowire is larger and accelerates the electron emission. 
Meanwhile, a higher transverse optical mode frequency 
results in a faster electron drift mobility. This makes 


the leakage current and polarization decay constant sig¬ 
nificantly thickness-dependent for thin wires. For thick 
nanowires, the soft mode frequency at a certain temper¬ 
ature approaches that of the bulk BaTiC >3 crystal. In 
this range, thickness affects the electric field, electron 
mobility, and leakage current little. It is worth mention¬ 
ing that the conclusion, that leakage current depends on 
thicknesses only for thin nanowires or films, is consistent 
with previous studies f6ll-[63j. 

In this study, the surface adsorbate OH plays the role 
of top electrode, which is the source of nonequilibrium 
charge carriers. This modified Schottky model, simu¬ 
lating leakage current in ferroelectric oxides, differs from 
the traditional Schottky emission model in many aspects. 
The electrons are localized in a 2 p orbital, and their 
wave vector distribution is now not based on Fernh-Dirac 
statistics. Different from the traditional Schottky emis¬ 
sion model, in which the thermal population of electrons 
leads to the temperature dependence of leakage current, 
temperature dependent leakage current is attributed to 
change of electric field through the nanowire and elec¬ 
tronic mobility accompanied with the hardening (or soft¬ 
ening) of the transverse optical mode in this model. 


V. CONCLUSION 

In summary, the depolarization process of BaTiC >3 
nanowires has been studied by both experiment and 
first principles calculation. We investigated the mech¬ 
anisms which govern the polarization decay and drew 
several principles. A new proposed theoretical model, 
which combines molecular orbital theory, quantum tun¬ 
neling and the modified Schottky equation, could explain 
successfully the general trends in the temperature and 
nanowire thickness dependent decay rates. Our study 
demonstrates that the surface adsorbate plays a signifi¬ 
cant role in stabilizing ferroelectricity and that depolar¬ 
ization is a process of charge leaking from the hydroxyl 
surface adsorbate to the gold substrate. 
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Appendix A: Derivation of equation (4) 


Appendix B: Derivation of equation (8) 

The derivation follows the idea in Ref. [H|. We con¬ 
sider that case that an electron with the wave vector k 
penetrates a BaTiC >3 (BTO) nanowire. 
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FIG. 5. Energy barrier diagram. 

As shown in FIG.0 The length of the penetration path 
is d. The energy of the electron is E (k) and the energy of 
the BaTiC >3 conduction band is V 2 . In the region of OH, 
the wave function is the plane wave including incident 
part and reflection part 

<pi (ar) = Ae ikx + Be~ ikx (x < 0) (Bl) 
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In the BaTiOa region, the Schrodinger equation is 


h 2 d 2 ip (x) 
2mo dx 2 


+ V (x) 
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(B3) 


When we consider the electrons around the conduction 
band, we can write the potential energy term into effec¬ 
tive mass 
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E BTO = V 2 ~E{k) (B5) 

where m* is the effective mass. The solution of the above 
Schrodinger equation is 

ip 2 (x) = Ce kx + De~ kx (0 < x < d) (B 6 ) 


kCe k ' d - k De~ kd = ikFe ikd (B13) 


From equation (B12) and (B13) we have, 

(V + ik) Fe ikd 


c = v ', 

2k'e k d 

(B14) 

(k — ik] Fe lkd 


D = V 2 , 

2k' e k d 

(B15) 

From (BIO)xifc+(Bll), we have 
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Substitute (B14) and (B15) into (B16) gives 
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With the analysis above, wave function in each region 
is summarized as 


The tunneling probability is given by 
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